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ABSTRACT 

In June-September 2006 the Be/X-ray binary EXO 2030+375 experienced the second giant outburst since its discovery. The source 
was shown to have a complicated pulse-averaged X-ray spectral continuum with possible evidence of cyclotron absorption features. 
In this paper we present the first pulse-phase resolved analysis of the broad band X-ray spectra of EXO 2030+375 obtained with the 
INTEGRAL observatory close to the maximum and during the decay phase of the giant outburst. We report a strong variability of 
the spectrum with pulse phase. Alternative spectral continuum models are discussed. The dependence of the spectral parameters on 
pulse phase during the maximum of the outburst and the evolution of the pulse profiles with time are qualitatively consistent with the 
pulsar's emission diagram changing from the fan-beam geometry close to the maximum of the outburst to a combination of pencil 
and fan beams (of comparable intesities) at the end of the decay phase. Evidence of a cyclotron absorption line around 63 keV at the 
pulse phase interval preceeding the main peak of the pulse profile is present in the spectrum obtained close to the maximum of the 
outburst. 
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1. Introduction 

■ The transient accreting pulsar EXO 2030+375 belongs to the 
' most common type of X-ray pulsar systems - the Be/X-ray bi- 
naries. Such systems form a subclass of high mass X-ray bina- 
ries. They consist of a pulsar and a Be (or Oe) companion, a 

^ main-sequence star of s pectral type B (o r O) that shows Balmer 
emission lines (see e.g. ISlettebaklll988 l. for a review). The line 
emission is believed to be associated with an equatorial outflow 
of material expelled from the rapidly rotating Be star that proba- 

■ bly forms a quasi-Keplerian d isk around its equator ( Hanuschik 
' Il996l:l0uirrenbach et al.lll997h . If the disk reaches a radius com- 
' parable to the periastron separation, then disk gas accreted by 

the neutron star can power a significant (and usually transient) 
X-ray source. 

Be/X-ray binaries typically show two types of outburst be- 
havior: 

1 . Normal (or type I) outbursts. They are characterized by rel- 
atively low X-ray luminosities Lx ~ 10^^-10^^ erg s"', low 
spin-up rates (if any), and recurrence at the orbital period (or 
its multiples). Such outbursts last from days to weeks and are 
associated with the periastron passages of the neutron star 

2. Giant (or type II) outbursts. They are characterized by higher 
X-ray luminosities Lx > 10^^ erg s"' and high spin-up rates. 
Such outbursts occur irregularly. They last several weeks and 
are not correlated with any particular orbital phase. The typ- 
ical time between outbursts is around several years. Giant 
outbursts are thought to stem from a dramatic expansion of 
the disk surrounding the Be star, leading to the formation of 
an accretion disk around the compact object. 

EXO 2030+375 is one of the best-studied Be/X-ray binaries. 
It was discove red with the EXOSAT satellite during a giant out- 
burst in 1985 dParmar et alJll989bl) . The optical companion of 



the pulsar is a BO Ve star identified by optical and i nfrared ob- 
serva t ions of the EXOSAT error circle (Motch & Jan ot-Pachecol 
Il987t iJanot-Pacheco et alj|1988t ICoe et al.iil98 8). The orbital 
period and eccentricity of the system are ~46 d and ~0.42, re- 
spectively (IWilson et al.ll2002i) . The period of X-ray pulsations 
is ~42 s. There were two giant outbursts in the history of obser- 
vations of EXO 2030+375. During the first one in 1985 (when 
the pulsar was discovered), the X-ray luminosity of the source 
reached a value of Li_2o keV ~ 2 x lO "^*^ ergs'' (assuming a 
distance of 7.1 kpc, I Wilson et al.l |2002|) . The spin frequency 
of the pulsar changed dramatically, with a spin-up timescale 
-P/P ~ 30 yr indicating the formation of an accretion disk 
around the neutron star. The second giant outburst took place in 
June-September 2006 dCorbet & Levind 120061: iKlochkov et"aLl 
SoOTI) and was again accompanied by a strong spin-up of the 
neutron star. The X-ray luminosity at the maximum of the out- 
burst was slightly lower tha n during the 1985 gian t outburst: 
ii-20kev ~ 1.2 X lO^^ergs-i dWilson & FingeJl2()06i) . 

During the 2006 giant outburst the source was observed sev- 
eral times with the INTEGRAL satellite. A preliminary pulse- 
averaged spectral analysis of som e of these observations are pre- 
sented in iKlochkov et al.l (l2007l) . A detailed analysis of pulse- 
averaged R XTE spectra o btaine d during the outburst was per- 
formed by IWilson et all (|2008). It has been shown that the 
pulse averaged X-ray continuum of the source has a com- 
plicated shape a nd cannot be mod eled by a simple power 
law/cutoff model. Wilson et a T'('2008') include an absorption line 
at -lOkeV (which they interp reted as a cyclotron line) in their 
spectral model, whereas [Klochkov et al.l (120071) have shown that 
the spectrum can be fitted equally well without the absorption 
line, but including a broad emission "bump" at ~15keV. 

Here we present for the first time pulse-phase resolved broad 
band (3-150keV) spectra of EXO 2030+375 during a giant out- 
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Table 1. Summary of observations. 
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Fig. 1. The RXTE/ASM light curve of EXO 2030+375 show- 
ing the 2006 giant outburst. Times of INTEGRAL observations 
analyzed in this work are marked with vertical lines. The obser- 
vations are referred to as Set 1, Set 2, and Set 3 in the text. 



burst. For our analysis we used all the available INTEGRAL 
data taken during the 2006 outburst. The X-ray continuum of 
the source shows strong variability with pulse phase, with some 
features present only at particular pulse-phase intervals. The de- 
scription of observations that we used is provided in Sect. |2] 
Details of data processing are described in Sect. [3] Sections |4] 
and |5] are devoted to the analysis of the data close to the maxi- 
mum and during the decay of the outburst, respectively. The re- 
sults are discussed in Sect. |6] and briefly summarized in Sect. |7] 



2. Observations 

The INTErn ational Gamma Ra y Astrophysics Laboratory 
(INTEGRAL. Wi nkler et"aD '2003') performed three pointed ob- 
servations of EXO 2030+375 during its giant outburst in June- 
September 2006. Two observations (on 19-20 August and 23- 
25 September) were done close to the maximum of the outburst. 
The third one (on 6-8 October) was performed at the end of the 
decay phase, when the X-ray luminosity dropped by a factor of 
~10 with respect to the maximum of the outburst. A part of the 
RXTE/ASM light curveQ including the giant outburst is shown 
in Fig. [U The times of the INTEGRAL observations are indi- 
cated and referred to as Set 1, Set 2, and Set 3 throughout the 
paper. 

Observations corresponding to Sets 1 and 2 were done 
when the source was at similar luminosity levels (Li-20kev ~ 
lO-'^ergs"') before and after the maximum of the outburst. 
A preliminary analysis of the two data sets was presented by 
iKlochkov et al.l (^007) who showed that spectral and timing 
characteristics of the source during the two observations were 
similar Due to the high X-ray flux, the data allowed us to 
perform a detailed pulse-phase resolved spectral analysis of 
EXO 2030+375 using twelve phase bins. During the observa- 
tions corresponding to Set 3, the X-ray luminosity of the source 
was at a level typical of normal (Type I) outbursts by this sys- 
tem (Li_20kev ~ 10^^ ergs"'). Spectra and pulse profiles of 
the source during these observations were considerably differ- 
ent with respect to those in Sets 1 and 2. In spite of a poorer 
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Mean 
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Table 2. The pulse periods and associated derivatives for three 
sets of observations analyzed in this work. 



Observation 


Reference 


P 


-dP/dt 




epoch (MJD) 


[s] 


[10-* s/s] 


Set 1 


53966.392169 


41.519270(5) 


2.758(20) 


Set 2 


53942.240864 


41.58084(3) 


2.96(21) 


Set 3 


54014.010200 


41.45425(6) 


0.22(23) 



Note: uncertainties in parentheses (68%) refer to the last digit(s). 



statistics, the data taken in Set 3 still allowed us to perform pulse- 
phase resolved spectral analysis, although with a coarser binning 
(four bins) in pulse phase. Table [1] contains the summary of the 
observations analyzed in this work. 



3. Data processing 

For our analysis we used t he data obtained wit h the instruments 
JB:S/fS GJ^I (20-300 key . fUbertini et al.l l2003h and JEM-X (3- 
30 keV, lLundetal.ll2003h onboard INTEGRAL. To perform the 
standard data reduction, the Off-line Science Analysis (OSA) 
software (version 6) was used (Cpurvoisier et al. 2003). To con- 
struct energy-resolved puls e profiles we used the software de- 
veloped at lASF, Palermo tMineo et al.ll2006l) . While extracting 
spectra and pulse profiles of EXO 2030+375, two other bright 
sources in the field of view of the INTEGRAL instruments, 
Cyg X-1 and Cyg X-3, were included to the extraction catalog 
used by the analysis software, which removes the contamination 
of EXO 2030+375. All other sources in the field of view are 
much weaker than EXO 2030+375 and therefore do not provide 
any noticeable contamination. The spectral analysis of the obser- 
vations was performed using the XSPEC v. 11. 3. 21 spectral fitting 
package dArnaudl 19961) . Following the OSA User Manual^ we 
added a systematic error to the final count rates at a level of 1 % 
for ISGRI and 2% for JEM-X to account for small-scale uncer- 
tainties in the response matrices of the respective instruments. 

To extract X-ray pulse profiles and to define pulse phase 
intervals for pulse phase resolved spectra, all photon arrival 
times were translated into the solar system barycenter and cor- 
rected for orbital m otion in the binary . The orbital parame- 
ters were taken from lWilson et alJ (l2005h : Poib = 46.0202(2) d, 
Tperi = MJD 5 1099 .43(2), a sin; = 238(2) It-sec, e = 0.416(1), 
cj = 210°. 8(4). The pulse periods and associated derivatives 
were determined for each of the three observational sets indi- 
vidually by employing initial epoch -folding and a subseq uent 
phase-conne c tion a nalysis similar to iFerrigno et al.l (l2006l) and 
iDeeter et al.l (Il981h . To perform phase-connection, we used a 
pattern-matching technique applied to well-defined pulse pro- 
files from a sufficient number of pulses. Any variation in the 
pulse shape inside each observation is marginal and does not 



' We used the results provided by the ASM/RXTE team. 
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affect our method. The determined periods and period deriva- 
tives are summarized in Table |2l Values corresponding to Sets 1 
and 2 are sligh t ly irn proved with respect to those reported in 
lKlochkovetal](l2007h . 



4. Maximum of the outburst 

Observations corresponding to Sets 1 and 2 were performed 
close to the maximum of the outburst (Fig. [T]i. During Set 1, 
which is roughly half as long as Set 2 (see Table [1]), the 
main target of observations was Cyg X-3 located 4° away from 
EXO 2030+375. Therefore, the statistics of the EXO 2030+375 
data obtained during Set 1 is noticeably worse than during Set 2. 
As a result, our pulse-phase resolved analysis of the maximum 
of the outburst is mainly driven by the data taken in Set 2. 

To explore spectral changes with pulse phase, we constructed 
pulse-phase resolved spectra and energy-resolved pulse profiles. 
While the former allow detailed study of spectral parameters as 
a function of pulse phase, the latter can be used to explore the 
intensity of different pulse-profile components (peaks, dips, etc.) 
in different energy ranges without relying on any particular spec- 
tral function. 

4. 1. Energy-resolved pulse profiles 

Figure |2] presents energy -resolved background-subtracted pulse 
profiles of the source obtained during the observations corre- 
sponding to Set 2. The corresponding pulse period and its deriva- 
tive, as well as the zero epoch, are provided in Table|2l The pulse 
profiles in the energy bands below 20keV were obtained from 
the JEM-X data, while the IBIS/ISGRl data were used for ener- 
gies above 20keV. Several components can be distinguished in 
the profiles: the main peak at phase ~0.2 followed by a fainter 
component at phase ~0.5 (which we will refer to as the "trail- 
ing shoulder") and two smaller peaks in the interpulse interval 
0.6-1 .0. The sharpest detail of the profile is a narrow dip preced- 
ing the main peak. Phase 0.0 was arbitrarily chosen to coincide 
with the dip in the average ISGRI pulse profile. As one can see, 
the shape of the profile changes smoothly with energy indicat- 
ing variations in the X-ray spectrum with pulse phase. The pulse 
fraction of the source determined as (Fmax - ■F,ni„)/(Fmax + ^'min) 
(where Fj^ax and F„^i„ are fluxes in the maximum and mini- 
mum of the pulse profile, respectively) as a function of energy is 
shown in Fig. [3] One can see a bump-like structure around ~10- 
20keV that is probably related to the continuum feature in this 
energy range (the "bump" or the absorption line, see Sects. [T]and 
14.2b . Above ~20keV, the pulse fraction seeply increases. 

4.2. Pulse-resolved spectra 

As one can conclude from the dependence of the pulse profile 
shape on energy, the X-ray spectrum of EXO 2030+375 clearly 
varies with pulse phase. We have therefore performed a separate 
analysis of the spectra accumulated in different pulse phase in- 
tervals. Phase binning was chosen to provide similar statistics of 
spectra in each bin and to have better phase resolution around the 
main peak where the most rapid spectral changes are expected. 
The binning was chosen a priori based solely on the shape of 
the pulse profile. No further adjustments were made after the 
spectra had been obtained. Figure|4]shows unfolded pulse phase- 
resolved spectra of the source (shifted vertically with respect to 
each other to avoid overlaps). The pulse phase is defined in the 
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Fig. 4. Pulse-resolved broad band X-ray spectra of 
EXO 2030+375 obtained in Set 2. Phase bins are indi- 
cated (phase zero is the same as in Fig.|2|. Values in parenthesis 
indicate the multiplicative factor applied to the flux in each 
spectrum to avoid overlaps. The solid line represents the fit of 
the spectra using Model I (see text for the model details). 



same way as in the previous section. Variability of the spectral 
continuum is clearly seen. 

As pointed out in Sect. [1] the pulse-averaged spectral con- 
tinuum of the source is rather complicated and cannot be mod- 
eled by any of the simple spectral functions (a power law mod- 
ified at higher energies by an exponential cutoff), which are 
usually used to fit spe ctra of accreting pul sars. In an attempt 
to model the spectrum IWilson et al.l ( 12008 ) modified a power 
law/c utoff model by a Gau ssian absorption line at ~10keV, 
while iKlochkov etal] ( |2007|) included a broad Gaussian emis- 
sion component at ~15 keV instead. In this work we tried to use 
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Fig. 2. Energy-resolved pulse profiles of EXO 2030+375 obtained with JEM-X (below 20keV) and IBIS/ISGRI (above 20keV) in 
Set 2 (where the data have the best statistics). 



both models to fit the pulse-resolved spectra. The reduced in 
all phase bins is very similar for the two models. This does not 
allow one to ultimately choose between the two spectral func- 
tions. The two models are referred to as Model I (with a "bump") 
and Model II (with an absorption line at lOkeV) throughout the 
paper. 

Below we describe the spectral functions corresponding to 
the two models. The common part of Models I and II is the power 
law/cutoff continuum smoothed at the cutoff energy by a third- 
order polynomial: 

(E-^, if £ < £cutoif - A£ 

£-r.exp(-^^), if £ > £,„,„ff + A£ 
AE^ + BE^ + CE + D, if ^cutoff - AE < E < 
< -Ecutoff + AE, 

where K is the normalization coefficient, E the photon energy; 
r, ^cutoff, Efoid, and AE are model parameters. Numerical co- 
efficients A, B, C, and D are chosen to obey the condition of 
continuity for the function and its first derivative at the points 
(^cutoff - AE) and (^cutoff + AE). Model I includes a broad 
Gaussian emission component (a "bump") around 15 keV: 

iModell - icont + Abump CXp \ — 7, ) , (Z) 



^cnnt — K ■ 



2(rl 

bump 



where ^bump and crbump are the energy and width of the "bump", 
and ^Tbump is the numerical constant describing the intensity 



of the component. In Model II, the power law/cutoff contin- 
uum /com is modified by a multiplicative absorption line with 
a Gaussian optical depth profile: 



^Modelll - ^cont ' exp < -Tiine exp 



(E - -Eline ) 

2crl 
line 



2^ 



(3) 



where Ziiine, criine, and Tiine are the centroid energy, width, 
and the central depth of the line, respectively. Additionally, 
we multiplied the functions /Model i and /Model ii by the factor 
exp[-A^HO"bf(£^)] describing the low-energy absorption by cold 
matter in the line of sight. Here, crbf(£') is the photoabsorp- 
tion cro ss-section per hydrogen atom for matter of cosmic abun- 
dances (iBaluciriska-Church & McCammonI Il992h used in the 
phabs model of XSPEC, and A^h is the equivalent hydrogen col- 
umn density. We also added a Gaussian emission line to model 
the iron fluorescence line at ~6.4keV. The latter, however, im- 
proves the fit only slightly. For all pulse-phase resolved spectra 
the value of A^h was fixed to 2 x 10^^ cm"^ (the average of the 
best-fit values found in all phase bins). It was also found that the 
energy of the "bump" in Model I does not change significantly 
with pulse phase being close to 15keV. It was, therefore, fixed 
to this value for all phase bins. 

Figures|5]and|6]show the most important spectral parameters 
obtained using the two models as a function of pulse phase. The 
data correspond to the observations performed in Set 2, which 
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Fig. 5. Best-fit parameters of Model I (see text) as a function of 
pulse phase obtained by fitting the spectra from Set 2. 

The dotted curve shows the averaged ISGRI (20-120 keV) pulse 
profile scaled and shifted vertically to match the plot. Vertical error 
bars correspond to 90%-uncertainties. 



have the best statistics. Vertical error bars correspond to 90%- 
uncertainties. Note that the pulse phase range shown in the fig- 
ures is from 0.0 to 1.5 so that some data points are repeated (the 
same is true for Fig.fTOli. 

The same kind of analysis was performed using the obser- 
vations corresponding to Set 1 (slightly before the maximum of 
the outburst). Due to poorer statistics, we had to use a coarser 
binning compared to Set 2. The spectral parameters of the pulse- 
resolved spectra were less constrained in this case. However, the 
behavior of the parameters agrees with what is found in Set 2. 
We also notice that the best-fit parameters of the pulse-averaged 
spectra taken i n Sets 1 and 2 are also consistent with each other 
(see Table 3 in lKlochkov et al.ll2007h . 

Residuals left after fitting the spectrum by Models I and II 
from phase bin 0.06-0.12 obtained in Set 2 reveal an absorption 
feature at ~63keV (see Fig.|2l). Inclusion of a Gaussian absorp- 
tion line at this energy flattens the residuals. The best-fit energy 
of the line is 63.6+^^ keV for Model I and 63. 3+^ !^ keV for Model 
II. We checked the presence and the energy of the feature using 
other possible continuum models, such as Fermi-Dirac cutoff 
([Tanaka 1986) and the so-called nega tive and positive power- 
law times Exponential model (NPEX, iMakishima et alJI 19991) . 
including a "bump" or an absorption line at lOkeV to match the 
continuum. It was found that the presence and the energy of the 
feature are independent of the choice of the spectral function. 
The F-test probability that the line is due to statis tical fluctua- 
tions is ~3xlO"^ (see however fProtassov et alj|2002l about non- 
applicability of the F-test to line-like features). 
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Fig. 6. Best-fit parameters of Model II (see text) as a function 
of pulse phase obtained by fitting the spectra from Set 2. As in 
Fig.|5] the dotted curve shows the ISGRI pulse profile. Vertical 
error bars correspond to 90%-uncertainties. 
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Fig. 7. Residuals of the spectrum from the phase bin 0.06-0.12 
fitted with Model I (top) and Model II (middle). The bottom 
panel shows the residuals after fitting the spectrum with Model 
I where a Gaussian absorption line at ~63.6keV is included (the 
corresponding fit with Model II looks very similar). 

5. Decay of the outburst 

5. 1 . Energy-resolved pulse profiles 

As mentioned in Sect. |2l during Set 3, INTEGRAL caught the 
source at the end of the outburst's decay phase when the X-ray 
luminosity was ~10 times lower than the maximum of the out- 
burst (Fig [Til. Due to lower statistics we used much coarser bin- 
ning in energy (for pulse profiles) and pulse phase (for pulse- 
phase resolved spectra) than in Set 2 (4 instead of 12 in both 
cases). The resulting pulse profiles are shown in Fig. [8] As 
before, the corresponding pulse period, its derivative, and the 
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Fig. 8. Energy-resolved pulse profiles of EXO 2030+375 ob- 
tained with JEM-X (below 20keV) and IBIS/ISGRI (above 
20keV) in Set 3 (at the end of the outburst). 



epoch zero can be found in Table |2] To choose the time of phase 
zero in accordance with the one used for Set 2 we used the de- 
pendence of the pulse profiles shape on the X-ray luminosity 
presented in Fig. 1 of Parmar et al.l JTQSQab . The dependence al- 
lows one to identify features in the profiles (e. g. the dip around 
pulse phase 0.6) in both observations (Sets 1 and 2) and choose 
the zero epoch for both observations consistently. 

The dependence of the pulse profile on energy in Set 3 
(Fig. ISJ is mainly characterized by a decreasing relative ampli- 
tude of the peak at phase ~0.9 with respect to the main peak 
(phase ~0.2) with energy. The sharp dip seen in the 2-lOkeV 
at phase ~0.6 almost disappears at higher energies. Generally, 
one can see that the shape of the pulse profiles is substantially 
different as compared to the profiles obtained close to the maxi- 
mum of the outburst at similar energies. In Sect. |6]we show that 
the changes could be qualitatively explained if assuming that a 
pencil beam component appears at the end of the outburst. 



5.2. Pulse-phase resolved spectra 

To produce pulse-phase resolved spectra from the data obtained 
in Set 3, we used four phase bins. They cover each of the two 
peaks and two phase intervals of the interpulse. The spectral 
continuum can be modeled well by the power law/cutoff func- 
tion provided by Eq. [T] without inclusion of the "bump" or the 
lOkeV absorption line. Like in Sect. 14.21 we modified the model 
at lower energies by photoabsorption (with A^h ~ 2 x 10^^ cm"^) 
and added a Gaussian line to model the iron fluorescence emis- 
sion at ~6.4keV. Figure |9] represents the unfolded pulse-phase 
resolved spectra fitted with the described model and shifted ver- 
tically with respect to each other in order to avoid overlaps. (This 
was done multiplying the flux in each spectrum by a numerical 
factor indicated in parentheses.) Significant variation in the con- 
tinuum is clearly seen. The pulse phase is defined in the same 
way as in Fig. [8] 




Energy [keV] 

Fig. 9. Pulse-phase resolved broad band X-ray spectra of 
EXO 2030+375 obtained in Set 3. Phase bins are indicated 
(phase zero is the same as in Figs.|8]l. Values in parentheses indi- 
cate the multiplicative factor applied to the flux in each spectrum 
to avoid overlaps. The solid line represents the fit of the spectra 
using the power law/cutoff model (see text for the model details). 
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Fig. 10. Best-fit spectral parameters of EXO 2030+375 obtained 
in Set 3 as a function of pulse phase. The dotted curve shows the 
averaged ISGRI (20-120 keV) pulse profile. Vertical error bars 
correspond to 90%-uncertainties. 



Figure \TU\ shows the photon index F and the folding energy 
iifoid as a function of pulse phase during Set 3. Vertical error bars 
correspond to 90%-uncertainties. 



6. Interpretation and discussion 

6. 1 . Pulse-period derivatives 

The values of pulse-period derivatives measured close to the 
maximum and at the end of the outburst (Table|2]l basically con- 
firm the statement in Klochkov et al. ( 2007) that the spin-up rate 
of the pulsar is proportional to its X-r ay flux as expected i n 
the simple accretion torque theory (e. g. lPringle & Reeslll972h . 
The observed dependence is similar to what was measured by 
EXOSAT during the 1985 giant outburst of EXO 2030+375 
(Parmar et al. 1989b). However, with only a few points, it is not 
possible to explore the correlation using more compUcated ac- 
cretion torque models. 
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6.2. Variation in fhe spectral continuum witli pulse phase 

The INTEGRAL observations have shown that the spectral con- 
tinuum of EXO 2030-1-375 is strongly pulse-phase dependent. 
Observations of the source close to the maximum of the out- 
burst allowed us to perform a detailed study of this dependence. 
A harder main peak and softer interpulse region observed in 
the source is a c ommon property of accreting pulsars (see e.g. 
iTsYgankov et al.ll2007L end references therein). It is usually ex- 
plained by assuming that, during the peak, we mostly see the 
comptonized photons coming from a hot region close to the foot- 
step of the accretion column while in the interpulse softer ra- 
diation scattered by the upper parts of the column is observed 
dBasko & Sunyaevlll976h . A closer look at Fig. |2] however, re- 
veals a more complicated picture. Both the main peak and its 
left flank increase in relative amplitude with energy, leading to 
a slight shift of the profile maximum towards an earlier pulse 
phase. This shows that the spectral continuum changes asym- 
metrically with respect to the main peak. The variation in spec- 
tral parameters with pulse phase (Figs. |5] and |6]l demonstrates 
the corresponding complex pulse-phase dependence of the spec- 
trum. The photon index F reaches a minimum (the hardest spec- 
trum) in the middle of the main peak (top panel in Figs. |5] and 
|6]l, while the maximum of the exponential folding energy Efou 
is shifted with respect to the peak towards earlier pulse phases. 
This contrasts to e.g. HerX-1, where both the maximum of 
£foid and the minimum of F are coincident with the main pulse 
dKlochkov etal . 2008). 

A possible qualitative explanation for the observed asymme- 
try of the spectral variation with pulse phase in EXO 2030-1-375 
can be the following. Since the luminosity of the source close 
to the maximum of the outburst, -lO"*^ erg s"', is well above 
the local E ddington limit (so-called critical luminosity, Lc ~ 
10^'^ erg s"'. iBasko & SunvaevllT976t IStaubert et alll2007h . it is 
expected that accreted matter is decelerated by radiation pres- 
sure leading to the formation of an optically thick accretion col- 
umn. The bulk of the emission in this case is expected to oc- 
cur as a fan beam (see als o the modeling o f the source's pulse 
profiles performed by Par mar et al.l Il989al) . Therefore, during 
the maximum of the profile (main peak) the angle between the 
column axis and the observer's line of sight has the highest 
value (the column is seen from the side). The observer is look- 
ing almost along the beam, seeing the photons coming from a 
Compton scattering region with large optical depth. This leads 
to the observed hard power law during the peak. Before the 
peak, when fifoid reaches a maximum, the direction of the line 
of sight might be closer to that of the column axis and, thus, of 
the magnetic field lines. Due to the dependence of the scattering 
cross-section on the angle between the photon direction and the 
magnetic field lines (iFlarding & Daughertvlll99lb . one expects 
that the photons, whose direction in this case is closer to that of 
the field lines, have experienced less scatterings and, therefore, 
originated deeper inside the accretion column where the temper- 
ature is higher The X-ray spectrum of these photons is expected 
to have larger Efow reflecting higher electron temperature but a 
softer power-law inde x due to lower Compt on scattering opti- 
cal depth (see e.g. Rvb icki & LightmanI 19791) . as observed. This 
explanation, however, requires that in the latter case the angle 
between the observer's line of sight and the column axis was not 
too small. Otherwise, depending on the geometry of the accre- 
tion flow, the column density along the line of sight will be very 
high, resulting in a higher optical depth. 

The changing of the pulse profile towards the end of the out- 
burst seems to confirm this picture. In the pulse profiles corre- 



sponding to the outburst's decay (Fig.[8j, one can see a new peak 
that appears at the phase interval preceeding the main pulse, i.e. 
where, according to our view, the observer's line of sight is clos- 
est to the magnetic field lines. At this pulse p hase one expects to 
see a pencil beam if the luminosity decreases jReig & Coell998l : 
iParmar et alj|l989al) . The observed peak, therefore, may corre- 
spond to the pencil beam component of the emission diagram, 
whose intensity is comaprable to that of the fan beam compo- 
nent at lower luminosity. Thus, both the pulse phase variation 
of the spectral continuum during the maximum of the outburst 
and the evolution of the pulse profile with luminosity fit the de- 
scribed picture where the emission diagram changes from a fan 
beam geometry close to the maximum of the outburst to a com- 
bination of a fan and a pencil beam at the end of the decay phase. 

The behavior of the spectral parameters with pulse phase in 
Set 3 is less clear due to lower statistics (only four phase bins are 
used). As in the maximum of the outburst, during the main peak 
of the profile, the spectrum is characterized by a hard F and a 
relatively low /ifoid (Fig.fTOl). During the peak around phase 0.0 
(pencil beam), /ifoid is higher. As before, this can be explained 
by photons from the pencil beam moving along the magnetic 
field lines originating deeper in the accretion column where the 
electron temperature is higher. On the other hand, the emission 
region is believed t o move closer to the sta r surface as the lumi- 
nosity decreases (Basko & Sunvaevl lI976). This means that the 
column density of gas above the emission region will be higher, 
leading to greater optical depth for Compton scattering. This can 
explain the hard photon index at this pulse phase. 

6.3. "Bump" versus absorption line at 10keV 

As discussed in Sect. 14.21 to model the spectral continuum of 
EXO 2030-1-375 in the maximum of the outburst, one has to mod- 
ify the standard power law/cutoff model either by a "bump" at 
~15 keV or by an absorption line at ~10keV. Both models pro- 
vide equally good fits of the pulse-phase resolved spectra. If in- 
terpreted as a cyclotron absorption feature (Wilson et al. 20081) , 
the absorption line at 1 ke V is observed at rather unusual energy 
that is much smaller than the exponential cutoff energy iicutoff- 
This contrast s with all other kno w cyclotron line sources (see 
e.g. Fig. 9 in Coburn et al.l l2002'). We point out here that a fea- 
ture around 8-15 keV (a "wiggle" or a "bump") i s observed in 
many accreting pulsars (e.g., 4U 1907, Her X- 1 . ICoburn et al.l 
2002), including those that do not otherwise exhibit a cyclotron 
fine (e.g. GS 1843-hOO, Coburn 2001). Thus, the interpretation 
of the absorption line at lOkeV in EXO 2030-1-375 as the cy- 
clotron resonant scattering feature might be hasty. On the other 
hand, the nature of the feature in this and other sources is still 
unclear Most probably it arises from modeling the spectral con- 
tinuum with a simple empirical function. For a proper modeling 
of the spectrum in the considered energy range (including the 
feature) one would need a proper theoretical model accounting 
for all relevant processes at the site of the X-ray emission. 

6.4. Evidence of a cyclotron line at ~63keV 

The X-ray spectrum of EXO 2030-1-375 taken close to the max- 
imum of the outburst in the narrow pulse-phase interval (0.06- 
0.12) preceeding the main peak shows evidence of an absorp- 
tion line around ~63keV (Fig.lT). If interpreted as a fundamen- 
tal cyclotron line, the corresponding magnetic field strength is 
B ^ 5x 10i2G(£cyc ~ 11.6x(B/10'2G) keV), which is one of 
the highest values among accreting pulsars. However, evidence 
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of a cyclotron line a t ~36keV ha s previ ously been reported for 
EXO 2030+375 by iReig & Cod (Il999h during a normal out- 
burst, so the line at ~63keV might well be the first harmonic 
rather than the fundamental line. It is known that the relative 
strength of the fundamental line and harmonics might vary sig- 
nificantly, sometimes making the fundamental line more difficult 
to detect than the harmonic (a good example is A 0535+26, see 
e.g. .Kendziorra et al...l994.) . 

7. Summary and conclusions 

We used INTEGRAL observations to study the pulse-phase de- 
pendence of the broad band X-ray spectrum of EXO 2030+375 
close to the maximum and during the decay of its 2006 giant 
outburst. This is the first pulse-phase resolved spectral study of 
the source. In all observations, significant pulse phase variability 
of the X-ray continuum was observed. 

Alternative spectral continuum models are discussed. We ar- 
gue that the interpretation of the feature at ~10keV as a cy- 
clotron absorption line proposed previously is questionable. 

Pulse-phase dependencies of the continuum parameters 
close to the maximum of the outburst, as well as the evolution of 
the pulse profiles from the maximum to the end of the outburst, is 
qualitatively consistent with the picture where the pulsar's emis- 
sion diagram changes from the fan-beam configuration close to 
the maximum of the outburst to a combination of pencil and fan 
beams (whose amplitudes are comparable) at the end of the de- 
cay phase. 

Evidence of an absorption line at ~63 keV is found during 
the maximum of the outburst at a narrow phase interval preceed- 
ing the main peak of the pulse profile. This feature can be inter- 
preted as the first harmonic of the previously reported cyclotron 
line at ~36 keV. 
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